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I. INTRODUCTION 



Unstable particles (such as the vector mesons) have an important role to play in low mass 
dilepton production resulting from heavy ion collisions and this has been a subject of several inves- 
tigations. Yet subtleties, such as the implications of the generalization of the Breit-Wigner formula 
for nonzero temperature and density, e.g. the question of collisional broadening, the role of Bose 
enhancement, etc., the possibility of the kinematic opening (or closing) of decay channels due to 
environmental effects, the problem of double counting through resonant and direct contributions, 
are often given insufficient emphasis. The present study attempts to point out these features using 
the rho and omega mesons as illustrative examples. The difference between the two versions of 
the Vector Meson Dominance Model in the present context is also presented. Effects of non-zero 
' temperature and density, through vector meson masses and decay widths, on dilepton spectra are 

studied, for concreteness within the framework of a Walecka-type model, though most of the basic 
fvj ' issues highlighted apply to other scenarios as well. 
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Heavy ion collisions at high energies produce matter far above the ground state providing thereby a rich arena 
for the study of hot hadronic matter, possibilities of chiral symmetry restoration, transition to a quark gluon plasma 
etc Q. However, hadronic signals are generally unsuitable for the task of uncovering the information on the underlying 
occurrences, since the history of strongly interacting particles entail layers of complicated dynamics which mask 
the basic issues. As such, electromagnetic signals, as manifested through emitted photon and dilepton spectra, are 
relatively cleaner, since electromagnetic quanta couple but weakly to hadronic matter. Final spectra exhibit resonance 
structures, which, in the low mass region, include the rho and omega mesons. Consequently, details of their creation, 
propagation and decay in the medium (or outside) are of paramount importance for the analysis of the resultant 
spectra observed. The general framework for such an investigation has been provided by Weldon Q through a 
. ' beautiful and lucid exposition on the Breit-Wigner (BW) formula at non-zero temperature and density. 

Accordingly section II, devoted to the underlying principles, begins with a quick review of the main results of 
Weldon's paper, followed by a brief discussion of the basic ideas of vector meson dominance (VMD) || which enables 
a phenomenological introduction of the coupling of photons to hadrons and hence to lepton pairs. We close our 
presentation of basics through an outline of a Walecka-type model Q which is merely used, in the present context, 
to provide a setting wherein we can formulate the hadronic scenario allowing us to estimate the differences that can 
accrue if insufficient emphasis is placed on the subtleties. We go on in section III to present the necessary ingredients 
to evaluate the dilepton emission rate from vector meson decays and the pion annihilation process in hot and dense 
hadronic surroundings. The last section is devoted to a discussion of results and conclusions. 



II. FORMALISM 



Ila. Generalised Breit-Wigner formula for unstable particles in a thermal bath 

Different species of hadrons in thermalised matter exist with equilibrium distributions determined by temperature 
T, the chemical potential \i and the statistics obeyed by that species. Consider an unstable hadron (R) in such a 
heat bath. If the decay products themselves are hadrons then they thermalise in the bath and no distinctive decay 
characteristic can possibly be discerned out. Thus we are interested in hadrons (R) that decay in the heat bath (one 
assumes that the collision volume can be so described) and decay into leptons and photons (described by the state 
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vector |/), say) that escape without thermalization. With q the total four momenta of the non-hadronic final state 
|/) the resonance peak should appear in the invariant mass (M) plot for say, the number of lepton pair events versus 
q 2 = M 2 at M = mR. The mass m R of the resonance R is the mass in the heat bath, theoretical estimation of which 
will of course depend on the model of hadronic interactions adopted. The central result of Weldon's paper is the 
generalization of the BW formula: 

dN f - (2i i n M2 r «^* r ™/( M ) m 



d 4 xd 4 q v ' 4tt 4 (M 2 - m 2 R ) 2 + {m R T 



tot \ 



with dNf the number of lepton pair events, say in the space-time and four momentum element d A xd i q 1 J being the 
spin of the resonance, r a ii->R is the formation width, Ttot is the total width and r^f > ^(M) is the partial decay width 
for off-shell R (i.e. of mass M) to go into the non-hadronic state |/). While this result is deceptively similar to the 
usual BW formula, it must be realised that the thermal distribution is implicitly contained through the 'entrance' 
width r aU ^ R as 

FalWR = exp^OE-^il' (2) 

and also through the interpretation of r to t which for a bosonic resonance R (our present concern) is given by the loss 
minus the gain, 

Ttot = r^^ai] — r a ii->R, (3) 

which is actually the rate at which particles equilibrate and relax to chemical equilibration. Here it is important 
to emphasize that the width of the invariant mass plots is related to the thermal damping rate r to t and this result 
generalises collision broadening treated by Van Vleck and Weisskopf [|| in the context of molecular spectroscopy. 

Reverting back to the problem at hand it is helpful to have a rough qualitative picture of decays for which the 
collision broadened BW is applicable, assuming that the thermalised hadron fluid lasts for a time Ty (after which it 
freezes out). The amplitude for a hadron to survive at time t can be modeled, albeit noncovariantly in the frame of 
the medium (bath), by 

A(t) = exp(-« E*t-m* R T tot t/2E*) < t < r f 

exp(-iEt-m R T™t/2E) t } < t < oo (4) 

where E and m R denote the energy and mass of the unstable particle and the asterisks represent the same quantities 
as modified by the medium. Particles for which r™£ Tt >> 1, the relevant portion of its history is from the early 
period (0 < t < Tf) as it is damped out at later times, and hence medium modifications determine its properties. On 
the other hand for particles with r™£ 77 << 1 and 1 > T/T to t >> T /rtot ^ * s * ne second term that dominates and in 
such cases the thermal effects on the mass and width are negligible. 

Another noteworthy feature is the fact that the decay width T^ 1 ^ * (M) , to be evaluated for an off-shell R, occurs 
in the formula and not /(M = m R ) though of course the point M = m R gets weighted most heavily (at the 
peak) due to the occurrence of the 'Breit-Wigner' denominator. However, for broad resonances this aspect does lead 
to some discernible differences as shall be illustrated later. 

Lastly, it is necessary to re-emphasize that for a particle which is sufficiently short lived to decay within the 
medium, the width of the dilepton spectra is actually the rate at which it equilibrates (r to t = Pr— >all ~ Taii^R), 
involving in principle various processes in which R participates. However it may be observed that elastic scattering 
does not enter into r to t as it cancels, contributing equally as it does to r^^an and r a n^R. This is however, in 
contradiction to earlier observations made in this context [|J. It may be borne in mind that although the elastic 
scattering contributes to the kinetic equilibrium it has no direct effect on the chemical equilibrium of the system 
while of course such elastic processes are of importance in phenomena such as viscosity etc. Indeed elastic scattering 
changes the momentum of the colliding particles but the nature of the particles remains unaltered and hence this 
process does not contribute to the decay life time in the bath. 



lib. Vector Meson Dominance 



The second important element in our framework is to have a robust phenomenological description of the electro- 
magnetism of hadrons (in particular rho and omega mesons) as this shall enter into the decay channels of interest. 
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Such a setting provided, for instance, by the Vector Meson Dominance (VMD) model proposed by Nambu and de- 
veloped by Sakurai ||, which assumes that the photon interacts with physical hadrons through vector mesons. Thus 
the cross section for the process ir + ir~ — » l + l~ (where I = e or p) will involve the coupling of the photon to the pion 
which is expressed in terms of the pion form factor F 1l (q 2 ) occurring in the matrix element of the electromagnetic 
current between pion states ( Tr(pf)\j p \ir(p) ) = F 1T (q 2 )(pf — p) M where the four momentum transfer is q = pi — p. While 
the photon can couple directly to the pion through its electromagnetic current eJ p = ie(ir~d p Tr + — Tr + d p iT~), the 
photon can also couple to the pion through a vector meson, which in this case, must also be an isovector. This is 
taken to be the rho meson. Based on such notions Sakurai enunciated the VMD model which has two formulations 
often referred to as VMD1 and VMD2 Q|. The photon and isovector meson part of the effective Lagrangian in the 
first representation is 

Cvmvoo = -\f^F^ - + \rn 2 p ^ Pli - g p ^J^ - eA^J^ - ^-F^ PlMJ , (5) 

where, p^ u is the field tensor for the rho field constructed analogously to the electromagnetic field tensor F^ , g p7rn 
may be determined from the decay p — > irir, and g p from fits to the process e + e~ — > tt + tt~. Thus here we have a 
direct photon-matter coupling as well as a photon-rho coupling which vanishes at q 2 = due to the occurrence of 
the derivatives in the last term above. This leads, after provision is made for the finite width of the unstable rho 
(emanating ostensibly from the imaginary part of the pion loop in the rho self energy), to the following expression for 
the rho dominated pion form factor, 

F VMm {q2) = j _ q 2 W 1 (6) 

g p q 2 - mf, + vm p Y p 

Note that the charge normalization constraint F 7r (q 2 = 0) = 1 is automatically built in. Furthermore, it may be 
remarked that often a gauge-like argument is advanced to the effect that rho couples universally (with the same 
strength) to all hadrons. However, the experimental fits reveal that universality is not exact; and indeed one finds 

= 1 + e, (7) 

with e = 0.2 Sakurai also outlined an alternative formulation (VMD2), which though not as elegant as the first 
(having for instance a photon mass-term in the Lagrangian), enjoys considerable popularity. Here 

1 11 6U1 2 P 2 

CvMve = — A F^F pu - jp^pp, + ^my Pfi - g p ^P^ J„ VA* + ^ m P A ^> (8) 

4 4 .9p ^9p 



and accordingly 



Furthermore 



where 



iem 2 

p — 7 vertex = . (9) 

9p 



F VMD2 { 2 } = m l ggZEl, (10) 

q 2 -m 2 +im p T p (q 2 ) g p 



and in order to maintain the condition F^^q 2 = 0) = 1, it is necessary here to impose the universality condition 
viz. g pjr7r = g pi whereas in VMD1 the charge normalization constraint is automatically maintained. Insertion of 
a momentum dependent width for the unstable vector meson in the VMD2 form factor maintains the condition 
F^(q 2 = 0) = 1. The momentum dependence originates on the one hand from the condition of the lowest mass state 
into which the rho meson can decay namely into two pions (hence the Heaviside theta function), and on the other 
hand due to the p-wave decay (therefore, the third power of the three momentum). We prefer VMD1 for reasons to 
be given later, the relevant results for VMD2 are also discussed as this version is used by several authors. Moving on 
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to the isoscalar part of the electromagnetic interactions of hadrons this shall analogously be taken to be dominated 
by the isoscalar vector meson w. The relevant part of the effective Lagrangian density involving u is taken to be |j||{J 

2 

/-relevant erU uj , ,u a . 9uirp ou u aa B \ 

L u = Lo h An H enva/30 oj a pr -it. (iz) 

The coupling of the omega to the photon (coefficient of uj^A p above) has long been considered to be approximately 
one third of that for the rho to the photon (coefficient of p^ A p in Eq.(g)), which yields reasonable agreement with the 
ratio of the observed partial widths T(p — > e + e~)/T(uj — > e + e _ ). Furthermore this is also supported by a recent 
QCD based study flC|| . The coupling g U7lp may be determined by using this term to calculate the observed w — > 7r° 7 
decay through the use of the rho-photon coupling already introduced in Eq.(^). However, the zero-ranged u> — 3tt 
vertex can also be obtained from the Lagrangian |pT[ , 

= / w3w e^w"e« fc a ,, 7r i a a ir i ^7r fc) (13) 

the latin indices referring to isospin. 



lie. The Walecka model and vector mesons in hot and dense hadronic matter 



The third ingredient needed, in order to discuss the characteristics of the rho and omega mesons in hadronic matter 
at a finite temperature (T) and density (jib) is to have an underlying model. Temperatures in the range ~ 150 — 200 
MeV and/or baryon densities nja few times nuclear matter density are of relevance. As a result the study of hadronic 
interactions leading to changes in their masses and decay widths under such conditions assumes great significance. 
Various investigations have addressed this issue over the past several years. Hatsuda and collaborators [jl2| and 
Brown [jl3| have used the QCD sum rules at finite temperature and density to study the effective masses of the 
hadrons. Brown and Rho JL4| also argued that requiring chiral symmetry (in particular addressing the QCD trace 
anomaly) yields an approximate scaling relation between various effective hadron masses, which implies, that all 
hadronic masses decrease with temperature. The gauged linear sigma model [ ^5[ , however, shows the opposite trend, 
i.e. m* increasing with temperature. In the present study we choose a Walecka-type model to be delineated below 
to provide the medium effects on the vector meson essentially for the sake of illustration. Many of our remarks and 
subtleties shall apply in some aspects to other scenario as well. 

The relevant interaction Lagrangian in the Walecka model, which we have considered, and comprising of the 
iso-scalar sigma, the rho, the omega and the nucleon, is given by, 

^relevant = g^^N - 9pNN (n^tN - i^l-Na ^tN d^j ■ ? - 9uN N N^N^ . (14) 

It may be observed that the rho has been taken to couple both minimally and through the Pauli tensor coupling to 
the nucleon. We shall adopt the value m a — 450 MeV for the sigma mass, and the coupling constants shall be taken 
to be g^NN ~ 10, g P NN ~ 2.6, k p ~ 6.1, k w = and g a ~ 7.4, chosen so as to reproduce the saturation density and 
the binding energy per nucleon in nuclear matter Jl6| . 

In the mean field approximation where the sigma meson field operator is replaced by its (classical) ground state 
expectation value (4> a ) ^ the Lagrangian ( |l4| ) immediately yields a medium dependent reduction in the nucleon 
mass. This is calculated from the nucleon propagator modified by a tadpole with a nucleon loop as its head and 
its tail emerging from the propagating nucleon line. Here we also include in the Relativistic Hartree Approximation 
(RHA), the properly renormalized contribution to the baryon self energy from the Dirac sea as well. This leads [|I||l7|] 
to a substantial reduction in the nucleon mass (m* N ) in the medium p8[ . It is in this setting that we consider the 
vector mesons. 

To compute in-medium meson propagators one solves Dyson's equation by essentially summing an infinite geo- 
metric series whose common ratio is the lowest order proper polarization which comprises of the nucleon loop (with 
in-medium mass m* N ) containing both the particle-hole (Fermi sea) and nucleon- antinucleon (Dirac sea) contributions. 
The effective mass of the vector meson (my) in nuclear matter is obtained by finding the value of energy go going to 
the limit q > for which the imaginary part of the propagator attains its maximum or equivalently by solving the 
full dispersion relation 

ql - q 2 - ™<v + Renf (T) (ffl,,q) + Ren F ( g 2 ) = 0, (15) 
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in the limit q — > in the region where Imll = 0. Hl(t) are the longitudinal (transverse) components of the in-medium 
self energy for thermal nucleon loop. 11^ is the vacuum self energy of the vector meson with modified nucleon mass 
due to sigma tadpole diagram. The expressions for Hl(t) and n> are given in the appendix of Ref. Jig] ]. 

It has been shown earlier Jl8|,[l^] that the change in the rho mass due to rho pion interaction is negligibly small at 
non-zero temperature and zero baryon density. In a different model calculation, it has been shown by Klingl et al pc| ] 
that to leading order in density the shift in the rho mass is very small. Therefore the change in rho meson mass due 
to p — tt — tt interaction is neglected here. However, at non-zero density the in-medium modification in the spectral 
function of the rho meson was studied [pl|-p3| by including the medium effects in the p — tt — tt vertex and the pion 
propagator in the delta-hole model. Since in this work we restrict our calculations within the realm of mean field 
theory (MFT) for internal lines, i.e the internal nucleon loop in the rho and omega self energy are modified due to 
tadpole diagram only, the inclusion of vertex corrections and modification of the pion propagator due to delta-nucleon 
hole excitation will take us beyond MFT and hence are not considered here for the sake of self consistency. Moreover, 
we do not include the delta baryon in the present work, as we have simply adopted a particular model for the sake of 
illustration. 



III. DILEPTON PRODUCTION 



The dilepton production rate due to processes occurring in a thermalised hadronic environment is obtained 
by folding the in-medium cross-section with the thermal distribution of the participants. In this article we consider 
dilepton production from pion annihilation (tt + tt~ — ► e + e~) and the rho decay (p — * e + e~). The thermal production 
rate per unit four-volume for lepton pairs is related to the imaginary part of the one-particle irreducible photon self 
energy by HJ||, 



dR _ tt/"Imn My (g) 
dS ~ l2TT 4 q 2 (ePi° - 1) 



(16) 



where, = (qo,q) is the four momentum of the virtual photon and n M „(g) is the one particle irreducible photon 
self energy. In the low invariant mass region, the pion annihilation channel is known to be the dominant one. The 
invariant mass distribution of the lepton pair in the case of pion annihilation is given by, 



dR 
dM 



M 3 
2(2tt)' 



(I - 4to 2 /M 2 ) j M T dM T dya(qo,ci) exp(-Mrcosh y/T) 



(17) 



Here, q = Mt cosh y, |q| = yffif — M 2 ; a(qo, q) is the cross-section for the pion annihilation calculated using VMD1 
and VMD2 described by Eqs. ([|) and (||) with appropriate modifications at finite temperature and density and is 
given by, 



c(g ,q) 



4 7TCT 

3 M 2 



y/1 - ±m 2 jM 2 yjl- Amf/NP (1 + 2mj/M 2 ) \F«{q , q)| 5 



(18) 



The (ft and q dependence of the pion form factors in the case of VMD1 and VMD2 originate from the real and 
imaginary part of the rho self energy in the medium. In the rest frame of the p meson, Eq. ( [T^ ) reduces to the well 
known form 



dR 
dM 



a{M) 



M 4 T Ki(nM/T) (1 - 4m^/M 2 ), 



(19) 



where K\ is the modified Bcsscl function, and M is the invariant mass of the lepton pair. 

In the same way, the invariant mass distribution of lepton pairs from the decays of vector mesons propagating 
with a momentum q is obtained using Eq.(Q), as 



dR *2J I 1 f 

— — = - 5 M M T dM T dy exp(-M T cosh y) 
dM 2tt J 



(M 2 - m 2 v + Ren) 2 
In the rest frame of the vector meson this reduces to 



-^ 2 r t 2 ot 



-pvac 



(20) 
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my V tot /n 



{MZ-mfY + m$T tot 



*2t^2 < n V L V -> e+ e- > 



(21) 



where r to t is defined by Eq.(g), and F 
vector particles. 



vac 



. (M) is the partial width for the leptonic decay mode for the off-shell 



IV. RESULTS AND DISCUSSIONS 



Within the ambit of the hadronic model adopted by us, the effect of finite temperature (T) and density (ub) 
on the self energies of the vector mesons reveals that the mass of the rho meson (to*) decreases more rapidly with 
increasing T and/or ns than that of the omega (to* ) as shown in Fig (|l|). One observes a small difference between the 
longitudinal(L) and transverse(T) modes in case of the omega meson but in case of the rho this splitting is negligible 
(attributable to the smaller vector coupling constant). We have observed that the quantity g 2 — q 2 along the dispersion 



curve remains almost constant^ 



,*2 



which is defined as % at q = on the mass hyperboloid) . This means that a 



simple pole approximation of the rho and omega propagator at k = my is good enough for our calculations. The 
splitting between the transverse and longitudinal components of the self energy of vector mesons with both vector 
and tensor interactions can be shown to be (see also ref. p(|), 



n, 



n L = 



^9VNN 



l_g2(^l)2 



k 2 dk d(cos ( 



2M' J J 
u cos 2 9 — v cos 6 + w 



M* 



C + 8k q |k| |q| cos 6 - 4k 2 q 2 cos 2 6 



Id + Id 



(22) 



where, 



3<Z 2 k 2 



g fco|k||q| ; w — 2fc 2 q 2 + q 2 k 2 — g 2 k 2 & C = q 4 — q^k^ and /d(/d) is the fermi distribution 



for nucleons (antinucleons). The imaginary part of the rho self energy is related to the probability of its survival in a 
medium. For a rho meson propagating in a medium with energy uj this is given by 



T{u J ) = 9 J^W\ S )- 
48-7T uj 



2T 



iy(s)v / ^ JTr 



In- 



1 - exp[-f (uj + W(s)Vuj' 2 - s)} 



1 — exp[- 



-4(0, 



w(s)Vu7 I ~Ts)} 



(23) 



where s = q 2 = uj 2 — q 2 and W(s) = y/T- 4to 2 /s. In the limit |q| — > 
in-medium decay width (with the decaying particle at rest in the medium) 



the above expression reduces to the 



9 p 7T 7T 

48tt 



loW 3 (uj) 



1 



(24) 



It is interesting to contrast the rate given by Eq. ( p3[ ) with the one obtained purely from considerations of relativistic 
time dilation from equation (p4). In this case the decay rate of an unstable particle in a general frame is given by, 



^Lor 



p rest frame 
P 



(25) 



From Fig. (g) it is clear that the decay rate of a particle propagating in a medium with momentum k cannot be 
obtained from its rest frame value by Lorentz transformation only. This is because of the preferred frame of the 
medium. 

It is well known that the naive Walecka model does not possess chiral symmetry, therefore, it is rather difficult to 
predict anything reliable on the pion mass. On the other hand, in models with chiral symmetry such as for instance 
the Nambu-Jona-Lasinio model, and the linear sigma model with nucleons, it is known that the pion mass is almost 
unchanged as long as it is in the Nambu-Goldstone phase. This is simply a consequence of the Nambu-Goldstone 
theorem in a medium [p7j . So we have adopted the usual approach of keeping the pion mass constant. 

The most significant process which contributes to the broadening of the omega in the thermal bath is uj tt «-> tttt, 
due to which the depletion rate of the omega increases substantially. This issue has also been discussed by Haglin || 
where the elastic channels involving the omega were also included. However, due to reasons mentioned earlier, we have 
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only considered the above (inelastic) channel to study the in-medium depletion of omega. Due to this process alone 
the decay rate of the omega comes out to be ~ 80 MeV. We now focus on an interesting possibility which becomes 
realisable (when rn* > m* + Tn w ), in that the decay to —* pn, which is closed in free space (as m u < m p + m w ), 
becomes an open channel. Thus whereas the omega meson which in usual circumstances decays through the three 
particle channel u> — > 3tt, could, given the appropriate environment, decay by the two particle (p7r) mode. Under 
these conditions the in-medium width for u> — * p tt is readily found by applying the finite temperature cutting rules 
to yield 



9ujirp 



327T77J 3 ,m?. 



A 3/2 ( 



ml, ml, 



) l + f(E^ + f(E p ) 



(26) 



where X(x, y, z) = x 2 



2{xy + yz + zx), arises from the phase space considerations, while / is the Bosc- 



Einstein distribution for the pions and the rho mesons in equilibrium. The coupling constant <? W7rp ~ 2 appearing 
in the expression T^^p^ can be deduced from the observed decay ui — > 7r°7, using the vector dominance model 
of Sakurai || for the pj vertex, taking the process to occur through a virtual rho converting to the photon. The 
concomitant three body process, lu — > 3tt, is estimated from the phcnomcnological effective term in the Lagrangian 
shown in Eq. (|l3|), averting for this purpose the Gell-Mann-Sharp- Wagner || model where the decay proceeds through 
a virtual rho lu — ► [pir] — * mm, in order to avoid the possibility of double counting when the threshold for the 
two body decay is crossed. In view of the above discussion one must include the processes lu <-> 3ir, lu <-> pir 
and urn — > irir as these are the most important among many other possible processes which can contribute to the 
broadening of omega in the medium. Taking all these ramifications into account the resulting width as a function 
of temperature at zero baryon density and at normal nuclear density (n° B ) is depicted in Fig. ([}]). We have noticed 
that the two body channel (uj — > pit) opens up in the former case at a temperature ~ 190 MeV, while in the latter 
situation (normal nuclear densities) this channel remains open even at zero temperature. We observe a value of ~ 86 
MeV for at a temperature of 200 MeV and normal nuclear density. As a result the lifetime of the omega reduces to 
~ 2.3 fm/c which should be compared to that of rho (~ 2.1 fm/c) under the same condition. This is in contradiction 
to the commonly held notion that the omega is too long-lived p8[ to convey any information on the fire ball in heavy 
ion collisions. 

While pointing out a potential phenomenon interesting in itself, we go on to emphasize the need in general to 
consider the possibility of the opening (or closing) of channels due to the effects of finite temperature or non-zero 
chemical potential. Modification of both the mass and width of the thermal omega due to the inclusion of anomalous 
interactions has also been discussed by Pisarski Ea,Bg|. 
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FIG. 1. Transverse and longitudinal dispersion relations of rho and omega mesons. The solid and dashed curves pertains to 
the transverse and longitudinal modes respectively. 
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FIG. 2. Decay rate of rho meson as a function of three momentum. Solid and dashed lines correspond to the in-medium 
width [Eq. (p^) ] and the boosted width [Eq. (^)] respectively. 
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FIG. 3. In-medium decay rate of omega meson comprising the processes u> — > pTY, lo — > 37r and ljtt — > iyit at ns = n% 
(dashed line) and ns = (solid line) as a function of temperature. 

It has been emphasized earlier that both the forward and backward processes should contribute to the probability 
of propagation of an unstable particle in a medium. This is manifested through the phase space factors which appear 
in the evaluation of such a quantity (which is neglected by some authors e.g. ^|). In this respect we discuss the term 
[1 + f{Ev) + f{Ep)} = (1 + / (E )) (l+f(E n )) — f(E p )f(E„), which occurs in the expression for the in-medium u> <-> pn 
width given by Eq. (p6|). This appears naturally in a calculation based on finite temperature field theory |p0| . Its 
physical significance resides in Bose-Enhancement (BE), which implies that the decay rate would increase because of 
stimulated emission in a gas already containing the decay products in equilibrium. Indeed the significant effect of this 
feature on the dilepton spectra has been demonstrated in a previous calculation [Q. This enhancement mechanism, of 
course, operates also for the p — > irir, toir — > it it and u> — > 3ir decays and has been incorporated in our calculations 
of decay width. 
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FIG. 4. Invariant mass distribution of lepton pairs from p — > e + e~ . The almost overlapping solid and dashed lines correspond 
to transverse and longitudinal rho with non zero three momentum. Dot-dashed line indicates the yield in the rest frame of the 
rho. 

However, as mentioned earlier, the hadronic decay modes of rho and omega in the fireball are not very informative 
and it is through their leptonic decay modes that they become experimentally 'visible'. Therefore, it is more relevant 
to examine the dilepton emission rate from rho in the medium using Eq. (^3|) which incorporates the generalised 
in- medium BW formula. Here we focus our attention on r to t- It is re-emphasized that though elastic processes such 
as p 7T «-> p 7r are there, they do not contribute to r to t • Indeed it should be borne in mind that though clastic collisions 
contribute to kinetic equilibration they do not contribute to the approach to chemical equilibrium, as indicated by 
Ttot = r_R^aii — Taii^R- For the case of the rho meson the processes p <-> nn and pit <-> u> are considered for the 
evaluation of r to t for rho. From phase space considerations it is clear that the mode p — » it it contributes dominantly 
to r to t- The resulting dilepton spectra from the decay of rho meson is shown in Fig. (||) at T=200 MeV and normal 
nuclear matter density. The notable feature here is the large shift of the rho peak towards lower invariant mass. This 
is due to the huge reduction in its mass [m* ~ 430 MeV). The solid and dashed curve show the resulting invariant 
mass distribution of the lepton pair from the decay of a rho propagating with momentum q in the thermal bath. 
The dot-dashed curve indicates the results when the decay p — > e + e~ is considered in its rest frame. A broader 
distribution in this case originates from the larger width of the rho in its rest frame as depicted in Fig. (||). However, 
it has been observed that the effect of collisional broadening due to p tt «-> won the dilepton yield from rho decay is 
insignificant. 




. Solid and dashed lines 



correspond to the case when T^O^f is evaluated for off-shell and on-shell rho (with non-zero q) respectively. 

We pass on to the discussion of the effect of off-shellness of the broad rho resonance on its dilepton decay mode 
as exemplified by the occurrence of T™^ l+l _(M) in Eq. (@) and again in Eq. ( po|) as contrasted width r^^ ;+; _ (jn* p ). 
Of course, in the narrow resonance limit when the BW structure reduces to a delta function peaked at m*, this effect 
is irrelevant. In Fig. (||) the solid curve indicates the dilepton yield when r v ^; e+e _ is evaluated at M (off-shell). The 
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on-shell result (dashed curve) shows a marked difference away from the rho peak (at the rho peak, of course, they 
must coincide). The off-shellness in rY? c , _ is calculated in the framework of VMD1. It is relevant to remark here 

\ — h 

that the in- medium 7 — p vertex is taken as em* /g P mr in VMD2 in the work of Li et al |31[ , which makes the coupling 
weaker due to the reduction in rho mass. 
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FIG. 6. Invariant mass distribution of lepton pairs from tt + ty~ — » e + e~ for VMD1 and VMD2. 

We devote this paragraph to the the process ir + ir~ — > e + e - , which could proceed through a photon coupled 
directly to the charge of the pion (ignoring its structure) and for q 2 ^ it would begin to see its structure which is 
modeled here by the intermediary rho meson. This separation is clearly manifested in VMD1 as can be seen from 
Eq. (JsJ) where the former mechanism is expressed through the occurrence of unity and latter exhibited by the rho-pole 
term. In VMD2, however, this feature is not at all manifest. A comparison of VMD1 and VMD2 results vis-a-vis 
dilepton yield from pion annihilation at temperature 200 MeV with tib = and tib = 2n^ is shown in Fig. (^) (the 
momentum dependence of the rho decay width is neglected here). The yield in the two cases is similar near the rho 
peak because the form factors around q 2 = M 2 ~ m* 2 become quite similar. However away from the peak M > m* 
the dilepton yield in the case of VMD1 dominates. For the reasons explained here and also in the introduction we 
have used VMD1 to evaluate the dilepton yield from pion annihilation and vector meson decays. 

Finally in order to see in what way the different issues considered above affect the low mass dilepton spectra, 
we plot in Fig. (0) the invariant mass distribution of lepton pairs from p — > e + e~ and irir — > e + e~ . We have used 
VMD1 to evaluate the contribution from pion annihilation. 

lO" 4 
f* 10" 6 
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■a 

« ICT? 
10" 8 

0.2 0.4 0.6 0.8 1.0 

M (GeV) 

FIG. 7. Invariant mass distribution of lepton pairs from the reaction tvtt — > e + e~ and the decay p — > e + e~ . Solid (dashed) 
line corresponds to the case when energy momentum dependence of the rho self energy is taken into account (ignored). 

In this work we have investigated the in-medium effects on dilepton production from pion annihilation and from 
the decay of unstable particles such as the rho meson (taken for illustration) . Subtleties arising due to the presence 
of a thermal bath, generalisation of the BW formula in the medium, collisional broadening and possibilities of double 
counting have been discussed. The variation of effective masses and decay widths of nucleons and vector mesons at 
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non-zero temperature and baryon density have been calculated within the framework of Walecka model. The BE 
effect in the decay width and the reduction in the mass of the rho meson is found to affect the dilepton emission rate 
quite substantially in the low invariant mass region. Moreover, in the presence of nuclear matter at finite temperature 
a substantial number of omega mesons could decay inside the reaction volume and thus can act as a viable probe 
for hot hadronic matter formed in relativistic heavy ion collisions. A comparison of the two versions of the vector 
dominance model has also been presented. 

Detailed measurement of photoproduction of lepton pairs should provide invaluable insights into the formation, 
propagation and decay of vector mesons inside the nuclear medium. Changes in the rho ( and also omega) masses 
would reflect directly in the dilepton invariant mass spectrum due to the quantum interference between rho and omega 
mediated processes in the photoproduction of lepton pairs (CEBAF Experiment). CERES collaboration |}2| has also 
planned to upgrade their experiment to improve the mass resolution such that rho and omega may be disentangled. 
Various aspects of the subtleties mentioned above on the observables with the inclusion of space time dynamics are 
under study, and are important for careful analysis of such experiments. 

Acknowledgement: We are grateful to Tetsuo Hatsuda for useful discussions. 
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